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PREFACE 


This  investigation  was  conducted  by  the  U.  S.  Army  Engineer  Water¬ 
ways  Experiment  Station  (WES),  Hydraulics  Laboratory  (HL),  from  October 
1976  to  October  1977  under  the  direction  of  Messrs.  H.  B.  Simmons,  Chief 
of  the  HL,  and  J.  L.  Grace,  Jr.,  Chief  of  the  Hydraulic  Structures  Divi¬ 
sion.  This  effort  was  supported  by  a  work  unit  (CWIS  No.  31042)  enti¬ 
tled  "Methods  of  Enhancing  Water  Quality"  of  the  Civil  Works  General  In¬ 
vestigations,  Environmental  Quality  Research  Area  and  Reservoir  Water 
Quality  Research  program  sponsored  by  the  Office,  Chief  of  Engineers. 
This  report  was  published  under  the  Environmental  and  Water  Quality 
Studies  (EWQOS)  program.  Dr.  J.  L.  Mahloch,  Environmental  Laboratory, 
was  Program  Manager  of  EWQOS. 

The  tests  were  conducted  under  the  supervision  of  Mr.  D.  G.  Fon- 
tane,  former  Chief  of  the  Reservoir  Water  Quality  Branch,  and  Dr.  E.  C. 
Tsivoglou,  Consultant.  Mr.  Steven  C.  Wilhelms,  in  partial  fulfillment 
of  the  requirements  for  a  Master  of  Science  degree,  Civil  Engineering, 
Mississippi  State  University  (MSU),  conducted  the  tests  and  prepared 
this  report.  Assisting  in  the  testing  and  analysis  were  Ms.  Cheryl 
Stevens  of  Dr.  Tsivoglou's  staff  and  Messrs.  M.  S.  Dortch,  C.  H.  Tate, 

B.  Loftis,  D.  H.  Merritt,  and  M.  E.  Neumann,  all  of  the  Hydraulic  Struc¬ 
tures  Division,  HL.  Dr.  V.  L.  Zitta  of  MSU  provided  guidance  and  tech¬ 
nical  assistance  on  the  preparation  and  presentation  of  this  report. 

COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE,  were  Com¬ 
manders  and  Directors  of  WES  during  this  investigation  and  report  prepa¬ 
ration.  Mr.  Fred  R.  Brown  was  Technical  Director. 
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TRACER  MEASUREMENT  OF  REAERATION:  APPLICATION  TO  HYDRAULIC  MODELS 


Hydraulic  Model  Investigation 


PART  I :  INTRODUCTION 

Background 

1.  Man  must  have  a  clean,  high  quality  water  supply  for  personal 
consumption,  growing  food,  and  processing  raw  materials.  Two  major 
sources  of  water  are  rivers  and  lakes,  thus  the  maintenance  of  high 
quality  in  these  resources  is  of  paramount  importance.  Clean,  high 
quality  water  is  the  most  important  element  to  aquatic  life  also,  since 
many  kinds  of  life  cannot  survive  in  low  quality  or  polluted  water. 

2.  Oxygen  is  the  one  necessary  element  for  survival  of  nearly  all 
life  forms.  Low  concentrations  of  oxygen  in  the  aquatic  environment  is 
a  hazardous  condition  for  many  kinds  of  life  that  exist  there.  Low  dis¬ 
solved  oxygen  (DO)  in  rivers  and  lakes  is  caused  by  the  natural  biologi¬ 
cal  processes  of  assimilating  organic  material.  It  is  important  to  note 
that  streams,  because  of  their  movement,  are  well  mixed  and  reaerate 
with  oxygen  from  the  atmosphere. 

3.  Consider  a  lake  in  early  spring;  the  water  is  isothermal  with 
a  constant  DO  content  from  surface  to  bottom.  As  summer  approaches, 
solar  energy  stratifies  the  lake  into  a  warm  upper  region,  the  epilim- 
nion,  and  a  cold  lower  region,  the  hypolimnion.  Wind,  which  had  been 
mixing  the  entire  lake  under  isothermal  conditions  and  causing  oxygen  to 
be  distributed  throughout  the  water  body,  can  now  only  mix  the  epilim- 
nion  since  the  hypolimnion  is  cold  and  heavy  and  remains  in  the  lower 
levels  of  the  lake.  Organic  material  in  the  hypolimnion  decomposes 
using  oxygen  and  since  mixing  does  not  distribute  atmospheric  oxygen  to 
the  lower  levels  of  the  lake,  the  hypolimnion  becomes  low  in  DO  or  even 
anaerobic.  Once  this  occurs  it  will  usually  remain  until  the  late  fall 
when  cooling  surface  temperatures  allow  destratification  and  complete 
mixing  of  the  water  body. 


4.  Many  man-made  lakes  release  water  through  outlets  located  very 
deep  in  the  lake.  This  results  in  low  DO  in  the  stream  below  the  lake 
during  the  stratification  season  unless  the  release  water  is  reaerated 
sufficiently  to  increase  the  oxygen  content.  There  is  great  concern  for 
providing  adequate  gas  transfer  to  enhance  the  downstream  DO.  Thus  it 
is  extremely  important  to  understand  the  gas  transfer  characteristics  of 
the  structures  which  release  water  from  lakes,  since  it  may  be  possible 
to  build  structures  that  promote  gas  transfer. 

5.  There  are  several  numerical  and  experimental  methods  of  evalu¬ 
ating  oxygen  uptake  in  streams  and  they  are  discussed  in  the  next  chap¬ 
ter.  One  method  for  accurately  evaluating  stream  reaeration  capacity 
was  developed  by  the  Federal  Water  Pollution  Control  Administration.  A 

gaseous  tracer,*  which  was  the  basis  for  the  procedure,  was  first  tested 

2 

and  demonstrated  in  one  river,  and  later  applied  in  a  series  of  direct 

3 

field  measurements  of  stream  reaeration.  The  development  of  the  gas¬ 
eous  tracer  technique  marked  a  significant  advance  in  the  study  of  the 

reaeration  process.  The  technique  has  been  successfully  applied  to  nu- 
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merous  streams  and  is  becoming  accepted  as  one  of  the  best  methods 

available  for  determining  stream  or  river  reaeration  rates. 

6.  There  are  still  many  unanswered  questions  pertaining  to  reaera¬ 
tion  or  gas  transfer.  All  of  the  field  studies  mentioned  above  have 
been  on  streams  flowing  in  natural  channels,  with  none  in  manmade  struc¬ 
tures  such  as  stilling  basins  below  dams,  powerhouse  tailraces,  or  other 
reservoir  outlet  works  structures. 

7.  One  method  of  determining  the  gas  transfer  characteristics  of 
an  existing  structure  is  to  use  the  tracer  technique  and  test  the  pro¬ 
totype  structure.  However,  this  can  be  very  costly  and  is  not  always 
possible.  Modification  of  the  stilling  basins  or  energy  dissipators  to 
improve  gas  transfer  characteristics  of  existing  structures  is  usually 
very  expensive.  If  the  project  is  in  the  planning  stage,  then  a  proto¬ 
type  does  not  exist  and  accurate  means  of  predicting  gas  transfer  char¬ 
acteristics  of  proposed  structures  are  not  available.  Because  of  these 
problems  in  evaluating  the  effects  of  design  changes  in  hydraulic 
structures,  it  would  be  advantageous  to  test  a  hydraulic  model  of  the 
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structure  to  determine  its  gas  transfer  characteristics  and  evaluate  the 
effectiveness  of  alternate  designs  or  structural  modifications. 


Objective 

8.  The  objective  of  this  study  was  to  demonstrate  and  evaluate 
the  applicability  of  the  tracer  technique  as  a  tool  for  use  in  physical 
models  to  determine  gas  transfer  characteristics  of  hydraulic 
structures . 


Scope 


9.  Tests  conducted  in  a  physical  hydraulic  model  of  a  spillway 
and  stilling  basin  were  used  to  evaluate  the  sensitivity  of  the  tracer 
technique.  Tests  were  conducted  with  and  without  a  flip  lip  on  the 
spillway  to  determine  the  ability  of  the  technique  to  detect  the  effects 
of  design  modifications  on  gas  transfer  characteristics.  Since  repro¬ 
ducibility  of  results  was  of  prime  importance,  two  replicate  tests  were 
performed  for  such  design  (see  Table  3  in  Part  IV). 


PART  II:  RELEVANT  THEORY 

Physics  of  Reaeration 

10.  Reaeration  can  be  considered  as  two  processes  working  to¬ 
gether,  molecular  diffusion  and  dispersion  due  to  turbulent  mixing. 

These  are  characterized  by  the  capture  of  oxygen  molecules  from  the  at¬ 
mosphere  at  the  air/water  interface  and  subsequent  distribution  of  the 
oxygen  throughout  the  volume  of  the  water  body.  Molecular  diffusion  is 
caused  by  the  inherent  kinetic  energy  possessed  by  the  oxygen  molecules, 
whereas  dispersion  due  to  turbulent  mixing  results  from  the  application 
of  external  forces  on  the  volume  elements  of  water. 

11.  If  a  group  of  dissolved  molecules  could  be  placed  in  a  quies¬ 
cent  water  body  without  disturbing  the  water,  the  dissolved  molecules 

would  gradually  spread  out  through  the  water  and  eventually  achieve  a 

3 

uniform  distribution  throughout  the  entire  water  body  (Figure  1  ).  They 
would  accomplish  this  in  totally  still  water  because  of  their  own  kine¬ 
tic  energy  associated  with  their  surrounding  temperature.  The  average 
kinetic  energy,  KE  ,  can  be  quantified  as  3/2  kT  ,  where  k  is  the 

Boltzmann  constant  and  T  is  absolute  temperature.  In  terms  of  mass 

2 

and  velocity,  the  kinetic  energy  is  1/2  MV  with  V  being  an  average 
velocity.  The  movement  of  the  dissolved  molecules,  which  is  entirely 


Figure  1.  Diffusion  in 
quiescent  water 
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random,  causes  the  molecules  to  uniformly  distribute  throughout  the 
water  body. 

12.  Fick's^  first  law  quantifies  the  process  of  molecular 
diffusion: 


where 


J  =  -D 


dc 
m  ds 


(1) 


J  =  net  flux  of  molecules  across  any  plane  in  the  water  volume, 
gm 

2 

sec  -  cm 


2 

D  =  coefficient  of  molecular  diffusion,  cm  /sec 
m 

dc  3 

^7  =  concentration  gradient  across  the  plane,  gm/cm  /cm 

dc 

The  force  driving  molecular  diffusion  is  the  concentration  gradient 

where  dc  is  the  difference  in  molecule  concentration  on  the  two  sides 
of  the  plane  and  ds  is  the  infinitesimal  distance  from  one  side  to  the 
other. 

13.  Einstein7  developed  an  expression  for  evaluating  the  molecular 
diffusion  coefficient: 


D 

m 


_  RT 

“  N  f 
o 


(2) 


where 

2  2 

R  =  universal  gas  constant,  gm  -  cm  /mole  -  °K  -  sec 
T  =  absolute  temperature,  °K 

23 

=  Avogadro's  number,  6.023  x  10  per  mole 
f  =  friction  factor,  gm/sec 


The  friction  factor  f  is  related  to  the  ability  of  the  surrounding 

g 

medium  to  impede  the  movement  of  the  diffusing  molecules.  Stokes 
showed  the  friction  factor  to  be  directly  related  to  the  viscosity  of 
the  medium  and  the  radius  of  spherical  particles  falling  freely  thr*.*«£h 
the  medium.  Hence,  the  molecular  diffusion  coefficient  is  a  function  cl 
the  absolute  temperature,  viscosity  of  the  fluid,  and  the  size  of  the 
diffusing  particle. 
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14.  To  further  examine  the  molecular  diffusion  process  related  to 
gas  transfer,  consider  a  completely  quiescent  isothermal  water  body. 
Assume  the  water  body  is  initially  void  of  oxygen  and  the  only  source  of 
oxygen  molecules  is  the  overlying  atmosphere.  A  very  short  time  later, 
oxygen  molecules  move  into  the  surface  layer  of  water  because  of  their 
kinetic  energy.  The  molecules  move  in  random  directions,  thus,  some 
escape  back  into  the  atmosphere  while  others  move  deeper  into  the  fluid 
body.  Since  the  diffusion  downward  is  impeded  by  the  viscosity  of  the 
water,  the  molecules  tend  to  collect  rapidly  in  the  uppermost  layers. 

15.  The  net  rate  of  gas  molecule  entry  into  the  water  is  equal  to 

the  rate  of  entry  from  above  (constant,  because  the  overlying  atmosphere 

has  a  constant  oxygen  partial  pressure) ,  minus  the  rate  of  escape  back 

to  the  atmosphere  (proportional  to  the  dissolved  oxygen  partial  pressure 

9 

in  the  upper  layers  according  to  Henry's  law  ).  Hence  the  rapid  accumu¬ 
lation  of  gas  molecules  in  the  surface  layers  causes  the  net  rate  of 
entry  to  quickly  become  very  small.  According  to  Fick's  law,  the  mass 
transfer  rate  through  a  plane  is  proportional  to  the  concentration  grad¬ 
ient  across  that  plane.  In  this  case  the  concentration  gradient  becomes 
very  small  resulting  in  a  very  small  mass  transfer  rate  across  the  air/ 
water  interface.  Thus,  the  lower  water  volume  elements  are  "starved" 
for  oxygen. 

16.  Consider  that  same  body  of  water  with  no  dissolved  oxygen  and 
the  only  source  of  oxygen  molecules  being  the  overlying  atmosphere.  How¬ 
ever,  this  time,  the  water  is  being  agitated  causing  turbulent  mixing 
(Figure  2).  Turbulent  mixing  and  convection  cause  water  element  1  to 
come  in  contact  with  the  atmosphere.  The  contact  time  may  be  very  short 
but  because  the  concentration  gradient  from  the  atmosphere  to  the  volume 
element  is  large  the  rate  of  mass  transfer  is  large  (Fick's  first  law). 
Hence,  element  1  acquires  a  relatively  large  amount  of  oxygen.  The  ele¬ 
ment  is  then  swept  away  from  the  surface  and  comes  in  contact  with  vol¬ 
ume  element  2  which  has  no  dissolved  oxygen.  Because  the  concentration 
gradient  across  the  interface  of  elements  1  and  2  is  large,  the  rate  of 
mass  transfer  is  large  and  element  2  acquires  a  significant  amount  of 
oxygen.  All  the  volume  elements  in  the  water  body  multiply  this 
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Figure  2.  Diffusion  and 
dispersion,  effect  of 
turbulent  mixing 


process;  thus,  it  is  obvious  that  dispersion  caused  by  turbulent  mixing 
greatly  speeds  gas  transfer. 

17.  The  turbulent  mixing  process  and  dispersion  have  been  mathe¬ 
matically  described  by  the  following  equation^’*** 


9c  u9c  v9c  w9c  _  /  92c  92c  92c  \ 

9t  9x  9y  9z  '^^2  9y2  ^2  j 


(3) 


where 

c  =  substance  concentration,  gm/cm^ 
t  =  time,  sec 

u  =  velocity  in  x-direction,  cm/sec 
v  =  velocity  in  y-direction,  cm/sec 
w  =  velocity  in  z-direction,  cm/sec 

2 

D  =  molecular  diffusion  coefficient,  cm  /sec 
m  ' 

In  turbulent  flow  the  velocities  and  concentration  are  expressed  as  a 
mean  u  value  plus  a  randomly  fluctuating  component  u'  (Figure  3). 

18.  For  example,  the  velocity  in  the  x-direction  is  u  =  u  +  u*  . 
Velocities  in  the  y-  and  z-directions  and  concentration  are  expressed 
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similarly.  Making  these  substitutions  in  Equation  3  and  taking  a  time 
average  of  each  term  the  following  equation  is  obtained: 


9c  u9c  v9c  w9c  SCu'c1 )  9(v* c’ )  9(w' c' ) 

9t  9x  9y  9z  9x  9y  9z 


(4) 


where  the  horizontal  bar,  i.e.  u,  9(u'c')  ,  denotes  a  time  average. 
The  terms 


9(u' c' )  9(v’c’ )  9(w’c')  , 

9x  9y  9z 

are  the  turbulent  mass  flux  per  unit  area  and  are  similar  to  "Reynolds 
stresses" 


9(u')2  9(u'v')  9(u'w' ) 

9x  9y  9z 


(x-direction) 


(6) 


as  found  in  the  turbulent  momentum  transport  equations. 

19.  Boussinesq^  theorized  that  these  stresses  could  be  related 
linearly  to  the  gradient  of  average  velocity, 


(u'v') 


(7) 
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where 


2 

£  =  turbulent  exchange  coefficient,  cm  /sec 
Similarly,  and  by  analogy**  with  Fick’s  first  law  (Equation  1),  those 
turbulent  dispersion  terms  have  been  related  to  the  gradient  of  average 
concentration, 

*  -Ex  H  (8> 

where 

2 

=  turbulent  dispersion  coefficient  in  x-direction,  cm  /sec 
20.  The  turbulent  transport  process  equation  can  be  written  as 

9c  u9c  v9c  w9c  _  _  9^c  9^c  9^c 

9t  9x  9y  9z  m  _  2  »  2  a  2 

3  9x  9y  9z 

*  s  (Ex  if )  *  h  (Ey  I)  *  h  (Ex  H)  w> 

In  most  turbulent  flow,  the  dispersion  coefficient,  E  ,  is  much  greater 
than  the  molecular  diffusion  coefficient,  D  .  Thus,  molecular  diffu- 
sion  typically  limits  the  rate  of  reaeration.  But  turbulent  dispersion 
and  convection  result  in  a  very  high  rate  of  molecular  diffusion  at  the 
surface  at  all  times.  Thus,  the  rate  of  reaeration  is  greatly  dependent 
on  turbulence  and  convection.  Turbulence,  in  reference  to  reaeration, 
has  a  special  meaning  related  to  the  rate  of  surface  replacement  and 
dispersion  of  water  volume  elements. 

Other  Theories  of  Gas  Transfer 


21.  The  preceding  theory  conflicts  with  the  "two  film"  theory  of 

gas  transfer  developed  by  Lewis  and  Whitman.  In  a  well-mixed  system  a 

stable  saturated  surface  layer  cannot  exist;  constant  surface  replace- 
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ment  precludes  this.  Danckwerts  "surface  renewal"  theory  is  more 
aligned  with  the  theory  discussed  in  previous  paragraphs.  According  to 
Danckwerts,  fluid  elements  are  exposed  to  the  atmosphere  for  varying  pe¬ 
riods  of  time  and  the  chance  of  an  element  being  swept  back  into  the 


r 


fluid  body  and  being  replaced  by  another  at  the  surface  is  independent 
of  its  "exposure  time."  He  theorized  that  the  rate  of  surface  area 
renewal  is  constant  for  a  given  system.  He  surmised  that  the  fraction 
of  surface  area  with  age  between  t  and  t  +  dt  is  equal  to  the  frac¬ 
tion  of  surface  area  which  previously  had  age  between  t  and  t  -  dt 
minus  the  portions  which  had  been  replaced  in  the  time  interval  dt  ; 
obtaining  the  distribution  of  fractional  area  to  be 

0Ct)  =  re"rt  (10) 

where 

6(t)  =  fractional  area  distribution  as  a  function  of  time,  per  sec 
r  =  surface  renewal  rate,  per  sec 

22.  According  to  this  theory,  the  surface  elements  absorb  oxygen 
molecules  up  to  the  instant  they  are  swept  away  from  the  surface.  Cou¬ 
pling  this  theory  with  Fick's  law,  the  average  absorption  rate  was  de¬ 
termined  to  be  a  function  of  the  surface  renewal  rate,  r  . 
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23.  Dobbins  developed  the  "film-penetration"  theory  of  gas 
transfer.  He  reasoned  that  the  concentration  gradient  within  the  fluid 
body  was  negligible  when  compared  to  the  gradient  in  the  surface  layer, 
hence  the  bulk  concentration  could  be  considered  constant  and  the  diffu¬ 
sion  equation  only  applied  at  the  surface.  Dobbins  adhered  to  the  con¬ 
cept  of  a  thin  film  whose  elements  absorbed  gas  by  molecular  diffusion 
before  being  replaced.  He  used  the  surface  age  distribution  obtained  by 
Danckwerts  and  determined  that  the  average  absorption  rate  of  oxygen 
molecules  was  a  function  of  the  surface  renewal  rate  and  the  thin  film 
thickness . 

24.  These  theories  allow  for  increased  gas  transfer  through  in¬ 
creased  surface  renewal  rates.  They  do  not  account  for  the  increased 
mixing  and  molecular  diffusion  between  water  elements  within  the  body  of 
water  upon  which  the  reaeration  rate  is  highly  dependent.  Further,  in 
most  complex  natural  or  nonlaboratory  systems,  it  is  very  difficult,  if 
not  impossible,  to  evaluate  the  surface  renewal  rate  and  the  thin  film 
thickness . 
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Predictive  Models  of  Reaeration 


25.  Numerous  models  have  been  developed  for  predicting  the  reaera 
tion  coefficient.  Those  models  may  be  classified  as  theory-based,  semi- 
empirical,  or  totally  empirical.  Theory-based  models  are  developed  from 
conceptual  models  of  gas  transfer.  Semiempirical  models  relate  the 
reaeration  coefficient  to  the  rate  of  energy  dissipation  or  longitudinal 
dispersion.  Totally  empirical  models  relate  the  reaeration  coefficient 
to  the  average  flow  velocity  and  average  depth  of  flow. 

26.  The  most  well  known  theory-based  prediction  equation  is  the 

14 

O' Connor-Dobbins  model.  It  was  obtained  from  the  derivation  of 
Dobbins’  "film  penetration"  theory  of  gas  transfer  and  relates  the  rate 
of  reaeration  to  the  liquid  film  diffusion  coefficient,  mean  forward 
velocity,  average  water  depth,  and  channel  slope  by  assuming  that  the 
surface  renewal  rate  was  equal  to  the  inverse  of  the  time  scale  of  the 
vertically  fluctuating  velocity  component, 


where 

v'  =  vertical  velocity  fluctuations,  cm/sec 
&  =  mixing  length,  cm 

Prandtls'  definition  of  mixing  length  and  a  logarithmic  velocity  distri¬ 
bution  provided  a  basis  for  evaluating  the  surface  renewal  rate  in  terms 
of  the  mean  depth  of  flow,  the  channel  slope,  and  the  gravitational  and 
Von-Karman  constants.  For  a  relatively  deep  channel  with  a  presumably 
small  velocity  gradient,  the  surface  renewal  rate  was  related  to  the 
average  forward  velocity  and  the  average  depth  of  flow.  The  form  of 
the  model  is 


va 

K,  a  \ 

2  Hb 


(12) 


where 


1<2  =  reaeration  rate  coefficient  for  the  stream,  per  sec 
V  =  mean  forward  velocity,  cm/sec 
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H  =  average  depth  of  flow,  cm 
a,b  =  experimental  constants 

O'Connor  and  Dobbins  determined  that  the  constants,  a  and  b  ,  were 
0.5  and  1.50,  respectively. 

27.  Many  of  the  empirical  prediction  equations  have  the  same  form 
as  the  O' Connor-Dobbins  model,  but  with  different  exponents  on  the  aver¬ 
age  velocity  and  average  depth  of  flow.  These  models,  such  as  Bennett- 

Rathbun,1'*  Churchill,  Elmore,  and  Buckingham,^  Langbein  and  Durum, ^ 
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and  Bansal,  were  developed  using  multiple  linear  regression  and 
several  hydraulic  parameters  of  open  channel  flow.  Average  velocity, 
average  width,  mean  depth,  and  channel  slope  are  some  of  these 
parameters.  Least  squares  criterion  was  the  usual  method  of  fitting  the 
regression  equation  to  the  investigator's  data. 

28.  Table  1  shows  empirical  model  exponents  indicating  similari¬ 
ties  and  differences  among  models  due  to  development  with  different  data 

sets.  The  linear  proportionality  constant  differs  for  each  model  also. 
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29.  Semiempirical  models  such  as  Krenkel-Orlob ,  Cadwallader- 

20  3 

McDonnell,  and  Tsivoglou-Wallace  relate  the  reaeration  rate  coef¬ 
ficient  to  the  average  forward  velocity,  V  ,  slope  of  the  energy 
gradient,  S  ,  and  mean  depth  of  flow,  H  .  The  form  of  these  models  is 


K2  a 


(VS)a 


(13) 


Table  2  shows  the  exponent  values  for  the  above-mentioned  models.  The 
exponents  differ  again  due  to  development  with  different  data  sets. 
Other  semiempirical  models  relating  the  reaeration  rate  to  average  ve¬ 
locity,  slope  of  the  energy  gradient,  and  average  depth  include 
Churchill,  Elmore,  and  Buckingham1^  and  Bennett-Rathbun. These  pre¬ 
dictive  models  have  the  following  form,  which  is  very  similar  to  the 
semiempirical  models  discussed  above: 


K2  a 


vasc 


(14) 
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Table  1 


Empirical  Prediction  Models 


V 

K,  or  \ 
2  Hb 


a 

b 

O'Connor-Dobbins  (theory-based) 

0.50 

1.50 

Bennett-Rathbun 

0.607 

1.689 

Churchill,  Elmore,  and  Buckingham 

0.969 

1.673 

Langbein  and  Durum 

1.00 

1.33 

Bansal 

0.60 

1.40 

Table  2 

Semiempirical  Prediction  Models 

a 

b 

Krenkel-Orlob 

0.408 

0.660 

Tsivoglou-Wallace 

1.000 

0 

Cadwallader-McDonnell 

0.500 

1.000 
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where  the  coefficients  a  ,  b  ,  and  c  are  2.695,  3.085,  and  -0.823, 
respectively,  for  the  Churchill  et  al.  equation  and  0.413,  1.408,  and 
0.273,  respectively,  for  the  Bennett-Rathbun  equation. 


Remarks 


30.  It  is  very  difficult  if  not  impossible  to  apply  the  theories 
of  gas  transfer  to  real  situations.  It  is  most  difficult  to  determine 
the  appropriate  values  for  the  coefficients  or  other  parameters  used  by 
these  concepts.  The  same  is  true  for  most  numerical  models.  As  men¬ 
tioned  for  the  O' Connor-Dobbins  model,  values  for  the  surface  renewal 

coefficient  and  the  thin  film  thickness  are  not  easily  ascertained. 

5  21  22 

Rathbun,  Zogorski  and  Faust,  and  Lau  have  presented  a  state- 

of-the-art  review  and  evaluation  of  most  techniques  used  to  determine 

reaeration  coefficients. 

31.  Some  of  the  predictive  models  use  average  velocity  and  depth; 
if  the  hydraulic  conditions  are  constant  for  the  stream  or  channel  reach, 
then  these  equations  can  be  used.  In  most  hydraulic  structures,  the 
hydraulic  conditions  change  radically  through  the  structure,  and  the 
concepts  of  average  velocity  and  average  depth  cannot  be  easily  applied. 
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PART  III:  METHODOLOGY 


Alternatives  to  the  Tracer  Technique 

32.  Upstream  and  downstream  DO  contents  can  be  used  to  compute 

the  reaeration  rate  coefficient  by  performing  a  DO  balance  as  developed 

23 

by  Streeter  and  Phelps.  This  is  accomplished  by  accounting  for  all 
the  DO  sources  and  sinks  in  the  stream  reach  except  atmospheric  reaera¬ 
tion.  Thus  the  difference  between  upstream  and  downstream  DO  contents 
(after  accounting  for  the  sources  and  sinks)  is  attributable  to 
reaeration. 

24  25 

33.  Another  method  is  the  disturbed-equilibrium  technique  ’ 
developed  at  the  Water  Pollution  Research  Laboratory.  This  technique 
uses  sodium  sulfite  and  a  catalyst  (cobalt  chloride)  to  chemically  de¬ 
oxygenate  the  water.  By  observing  the  DO  uptake  in  the  stream  reach  and 
assuming  that  other  DO  sources  and  sinks  are  constant  during  measure¬ 
ments,  the  reaeration  coefficient  can  be  computed.  DO  uptake  measure¬ 
ments  made  while  using  this  method  are  presented  herein. 

Tracer  Technique 

34.  The  gaseous  tracer  method  for  determining  reaeration  uses 
three  tracers  simultaneously.  Krypton-85  (kr-85),  as  a  dissolved  gas, 
is  the  tracer  for  DO;  tritium,  as  tritiated  water  molecules,  is  the 
tracer  for  dispersion;  rhodamine-WT  fluorescent  dye  provides  an  ac¬ 
curate  measure  of  the  time  of  travel  and  indicates  when  to  sample 

for  the  two  radioactive  tracers. 

35.  To  understand  the  tracer  method,  the  equivalence  of  two  pro¬ 
cesses  must  be  made  clear.  The  absorption  of  atmospheric  gases  (pri¬ 
marily  oxygen)  into  water  and  the  desorption  of  tracer  gas  from  water  to 
the  atmosphere  are  equivalent.  In  either  case,  the  driving  force  for 

gas  transfer  is  the  difference  between  the  partial  pressures  of  the  gas 

9 

in  the  atmosphere  and  in  the  water.  Henry's  law  states  that  a  gas  will 
attain  a  saturation  concentration  in  a  liquid  that  is  proportional  to 
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the  partial  pressure  of  the  gas  in  the  overlying  atmosphere.  Thus, 
there  is  a  force  causing  molecular  diffusion  of  gas  molecules  into  the 
liquid  until  this  saturation  concentration  is  attained  equalizing  the 
partial  pressures  in  the  liquid  and  the  atmosphere.  For  oxygen,  a  mea¬ 
sure  of  the  strength  of  the  force,  called  "saturation  deficit,"  D  ,  is 
the  difference  between  the  oxygen  concentration  at  the  point  of  satura¬ 
tion,  ,  and  the  concentration  of  oxygen,  C  ,  which  actually  exists 

in  the  water  or 

D  =  (C  -  C)  (15) 

s  ox 


As  long  as  the  concentration  of  oxygen  in  the  water  is  less  than  the 
saturation  concentration,  there  will  be  a  net  movement  of  oxygen  from 
the  atmosphere  to  the  water.  This  process  usually  recognized  as  a  first 
order  reaction  process  can  be  derived  from  Fick's  law  and  is  represented 
by 


dD 

dt 


-K  D 


ox 


(16) 


where  t  is  time  and  K  is  the  reaeration  rate  coefficient  for  flow- 

ox 

ing  water.  This  equation  simply  states  that  the  rate  of  change  of  the 
saturation  deficit  at  any  time  is  proportional  to  the  deficit  at  that 
time,  or  the  greater  the  deficit,  the  greater  the  reaeration.  The 
magnitude  of  the  proportionality  constant,  Kqx  ,  is  dependent  particu¬ 
larly  upon  the  intensity  of  turbulent  mixing  in  the  system. 

36.  If  there  were  no  factors  other  than  turbulent  mixing  affect¬ 
ing  the  oxygen  concentration,  integrating  Equation  16  would  provide  a 
means  for  determining  the  proportionality  constant,  Kqx  ,  through  the 
relationship 

-K  t 

D  =  D  e  0X  (17) 

o 


where 

D  =  saturation  deficit  at  some  later  time,  t 

D  =  saturation  deficit  at  some  initial  time  (t  =  o) 
o 

There  are  many  chemical  and  biological  processes  which  affect  DO  in 
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1 


natural  systems,  and  therein  lies  the  value  of  using  the  inert  gas 
tracer  method.  The  inert  gas  is  not  affected  by  chemical  or  biological 
processes . 

37.  The  desorption  process  of  krypton  from  the  water  is  equiva¬ 
lent  to  the  absorption  process  of  oxygen  from  the  atmosphere.  The  driv¬ 
ing  force  in  the  desorption  of  the  tracer  gas  is  the  difference  between 
the  partial  pressure  of  the  krypton  in  the  water  and  the  partial  pres¬ 
sure  of  the  krypton  in  the  atmosphere.  For  all  practical  purposes,  the 
partial  pressure  and  thus  the  concentration  of  krypton  in  the  atmosphere 
is  zero.  Therefore,  tracer  gas  in  the  water  will  be  steadily  lost  to 
the  atmosphere.  A  measure  of  the  strength  of  the  force  causing  gas  loss 
is  simply  the  concentration  of  tracer  gas  in  the  water.  This  process 
can  be  represented  by 


where 


C  = 


C  e 
o 


-k.  t 

kr 


(18) 


C  =  concentration  of  kr-85  remaining  in  the  water  at  some 
later  time,  t 

Cq  =  concentration  of  kr-85  at  some  initial  time  (t  =  0) 

=  gas  exchange  coefficient  for  krypton  in  flowing  water 
Since  krypton  gas  is  inert,  it  is  not  subject  to  the  chemical  and  bio¬ 
logical  processes  that  affect  oxygen.  This  fact  makes  it  possible  to 
compute,  through  Equation  18,  a  gas  exchange  coefficient  for  krypton, 

Kfcr  ,  which  reflects  only  the  turbulent  mixing  process  and  direct 
physical  gas  transfer. 

38.  It  has  been  shown  that  the  ratio  of  exchange  coefficients  for 
these  two  gases  is  equal  to  a  constant,1  that  is, 


where 


0.83 


(19) 


K^r  =  exchange  coefficient  for  krypton 
Kqx  =  exchange  coefficient  for  oxygen 
This  relationship  is  not  affected  by  temperature  (within  the  range  of 
interest),  degree  of  turbulent  mixing,  or  the  direction  of  gas  transfer. 


20 


This  makes  possible  the  calculation  of  from  Equation  18  and 

subsequent  determination  of  Kqx  with  Equation  19. 

39.  Consider  two  points  A  and  B  which  lie  on  a  natural  system 
such  as  a  stream  or  channel.  Let  A  be  the  upstream  point  and  let  a 
quantity  of  dissolved  krypton  gas  be  introduced  upstream  of  A  .  If  the 
tracer  gas  were  introduced  uniformly  across  the  cross-sectional  flow 
area  and  there  were  no  vertical  or  horizontal  velocity  gradients  in  the 
flow  causing  dispersion,  and  if  there  were  no  tributaries  to  cause  dilu¬ 
tion,  then  the  exchange  coefficient  of  krypton-85,  ,  between  points 

A  and  B  could  be  calculated  from  Equation  18  in  the  form 


where 


:B  "Kkrt 
r-  =  e 


(20) 


CA,Cg  =  krypton-85  gas  concentrations  at  A  and  B 
t  =  time  of  travel  from  A  to  B 

40.  However,  since  dilution  and  dispersion  are  present,  they  must 
be  considered.  Even  if  it  were  possible,  direct  measurement  of  disper¬ 
sion  and  dilution  is  not  necessary;  a  correction  may  be  applied  to  Equa¬ 
tion  20  by  using  the  second  radioactive  tracer,  tritium. 

41.  Tritium,  in  the  form  of  tritiated  water  molecules,  is  re¬ 
leased  simultaneously  with  the  krypton  gas.  The  tritium  concentration 
decreases  from  point  A  to  point  B  since  it  is  subjected  to  the  dis¬ 
persion  and  dilution  of  the  reach.  Thus,  it  provides  an  accurate  mea¬ 
sure  of  the  dispersion  and  dilution.  Since  the  tritium  is  in  the  form 
of  tritiated  water  molecules,  it  is  not  adsorbed  on  the  channel  bottom 
or  sides  or  otherwise  lost  in  any  significant  amount.  Because  the  trac¬ 
ers  are  released  simultaneously  the  krypton-85  undergoes  the  same  dis¬ 
persion  and  dilution  as  the  tritiated  water. 

42.  The  observed  concentrations  of  tritium  therefore  provide  an 

accurate  correction  for  the  effects  of  dispersion  and  dilution  which  can 
be  applied  to  the  krypton-85  gas  concentrations,  and  the  krypton  ex¬ 
change  coefficient,  ,  can  be  calculated  by 


21 


-  v 


f 


where 


(21) 


ratios  of  krypton  concentrations  to  tritum 
concentrations  at  points  A  and  B 


t  =  time  of  travel  from  A  to  B 


Decreases  in  tracer  concentrations  due  to  natural  decay  need  not  be 
considered  since  the  isotopes'  half-lives  are  usually  much  greater  than 
the  travel  time. 


Model  Facility 

43.  The  model  used  in  this  study  was  a  l:20-scale  simulation  of  a 
single  bay  of  the  Lower  Monumental  Dam  spillway  located  on  the  Snake 
River  in  Oregon.  The  model  was  designed  and  constructed  for  hydraulic 
similitude  based  on  Froudian  criteria.  It  reproduced  a  short  portion  of 
forebay  and  a  single  bay  of  the  spillway  section  of  the  dam  with  the 
tainter  gate  for  flow  control.  The  stilling  basin  and  a  portion  of  the 
exit  channel  were  included  in  the  model  (Figure  4). 

44.  The  two  spillway  configurations  tested  were  with  and  without 
a  flip  lip  on  the  spillway  face  (Figure  5).  The  flip  lip  deflected  flow 
along  and  near  the  surface  of  the  tailwater  (Figure  6)  instead  of  allow¬ 
ing  the  water  to  plunge  to  the  bottom  of  the  stilling  basin  (Figure  7). 
This  subjects  the  air-entrained  flow  to  less  pressure  and  greater  expo¬ 
sure  to  the  atmosphere  which  results  in  greater  opportunity  for  gas 
transfer.  The  flip  lip  was  located  at  el  132.28*  which  was  below  the 
tailwater  surface  for  the  flow  which  was  tested. 

*  All  elevations  (el)  cited  herein  are  in  metres  referred  to  mean  sea 
level . 
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Figure  6.  Flow  conditions  with  flip  lip 
(15,000  cfs  =  425  cms) 
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Figure  7.  Flow  conditions  without  flip  lip 
(15,000  cfs  =  425  cms) 
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45.  Water  was  supplied  to  the  model  by  a  51-cm-diameter  pipe  with 
a  gate  valve  and  measured  with  a  venturi  meter.  The  flow  rate  for  all 
tests  was  0.238  cms  in  the  model  which  is  equivalent  to  425  cms  in  one 
bay  of  the  prototype. 
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PART  IV:  TESTING 


Apparatus  and  Equipment 

46.  Because  of  the  short  flow  time  through  the  model  and  quick 

2  3 

passage  of  the  dye  cloud,  previous  sample  collection  methods  *  could 
not  be  used.  A  high  speed  sampling  system  was  developed  for  use  at  both 
the  upstream  and  downstream  sampling  stations  (Figure  8).  A  submerged 


Figure  8.  High  speed  sampling  system 

pump  supplied  water  to  a  manifold  which  routed  the  flow  to  15  sample 
bottles  and  a  continuous  flow  fluorometer  that  was  equipped  with  a 
recorder.  Flow  to  the  sample  bottles  and  fluorometer  was  controlled 
with  adjustable  pinchcocks  and  was  calibrated  to  a  rate  of  10  mJfc/sec. 
Electromechanical  solenoid  valves  in  the  flow  lines  from  the  manifold 
to  the  sample  bottles  were  used  to  shut  off  flow  to  the  sample  bottles 
at  a  prescribed  time.  Water  was  continuously  supplied  to  the  sample 
bottles  and  fluorometer.  Excess  water  to  the  sample  bottles  was 
continuously  spilled  (Figure  9).  As  the  dose  passed  the  sampling  sta¬ 
tion,  part  of  it  was  pumped  through  the  manifold  to  the  sample  bottles. 
The  solenoid  valves  were  closed,  one  at  a  time,  at  a  prescribed  interval, 


thereby  shutting  off  the  flow  to  each  bottle  and  capturing  a  discrete 
sample  in  time.  Figure  10  is  a  schematic  of  the  sampling  system.  The 
two  systems  used  in  the  upstream  and  downstream  approaches  were  con¬ 
structed  to  be  identical,  with  the  length  and  volume  of  all  water¬ 
carrying  tubing  and  parts  minimized. 

47.  The  flow  was  instantaneously  dosed  with  the  radioactive  trac¬ 
ers  and  fluorescent  dye  by  crushing  the  glass  dose  container  with  a 
"bottle-breaker."  The  dose  container  was  "instantly"  crushed  by  strik¬ 
ing  a  plunger  with  a  hammer. 

Procedure 


48.  At  the  beginning  of  each  test  the  model  was  checked  for  cor¬ 
rect  flow  and  water  surface  elevations  and  the  hydraulic  situation  al¬ 
lowed  to  stabilize.  The  sampling  systems  which  were  located  at  the  two 
sampling  stations  (Figure  11)  were  put  into  operation.  Once  the  sampling 
systems,  the  f luorometers ,  and  the  strip  chart  recorder  were  ready,  the 
flow  was  dosed  with  the  tracer  mixture  (14.0  millicuries  of  krypton-85, 
6.0  millicuries  of  tritium,  and  fluorescent  dye)  upstream  of  the  model. 
Shortly  after  the  leading  edge  of  the  dye  cloud  passed  the  upstream 


FROM  OTHER  FLUOROMETER 


dual  channel  strip  chart 


SUBMERSIBLE  PUMP 

Figure  10.  Sampling  system  schematic 


sampling  station  (indicated  by  the  fluorometer  and  recorder),  sampling 
was  initiated  and  marked  on  the  chart  of  the  recorder.  A  sample  was 
taken  every  2  sec.  Every  fifth  sample  was  indicated  on  the  chart  to 
help  assure  analysis  of  the  correct  samples.  The  same  procedure  was 
followed  as  the  dose  passed  the  downstream  sampling  station. 

49.  As  soon  as  the  tracer  cloud  had  passed  and  all  fifteen  sam¬ 
ples  had  been  collected,  the  bottle-filling  tubes  were  carefully 
withdrawn,  and  the  bottles  were  tightly  capped  and  taped  shut  with 


FRO* I  SUMP 


Figure  11.  Sampling  station  location 

electrical  tape.  The  capped  samples  were  submerged  in  water  kept  at  a  constant 
temperature  for  transport  to  the  laboratory  for  analysis. 

Analysis 

50.  As  has  been  stated,  the  fluorescent  dye  tracer  was  used  to 
determine  when  the  dose  passed  the  sampling  stations.  With  the  sampling 
technique  used  in  these  tests,  samples  were  taken  over  a  large  portion 
of  the  dye  cloud  assuring  the  collection  of  samples  around  the  time  of 
peak  dye  concentration  and  peak  radioactive  tracer  content. 

51.  The  travel  time  through  the  model  was  determined  from  the  re¬ 
corded  fluorescent  dye  curves.  From  center-of-mass  (upstream  curve)  to 
center-of-mass  (downstream  curve)  could  be  considered  as  the  travel  time 
between  the  sampling  stations.  However,  collecting  samples  for  analysis 
near  the  center-of-mass  of  the  dye  cloud  was  not  practical  since  the 
center-of-mass  was  indeterminable  until  after  the  dye  cloud  had  passed 
the  sampling  station.  Further,  the  krypton  gas  transfer  coefficient, 

K^r  ,  computed  with  Equation  21  is  only  valid  when  used  with  the  travel 


time  with  which  it  was  calculated.  By  examining  Equation  21  it  can  be 
seen  that  the  unknown  in  the  equation  is  actually  the  product  of  K^t  . 
Substituting  any  observed  time  of  travel  would  result  in  a  different  ex¬ 
change  coefficient,  but  the  product  would  remain  the  same. 

52.  For  analytical  purposes,  the  samples  of  principal  interest 
were  those  taken  when  the  tracers  were  at  their  peak  concentrations. 
These  samples  were  chosen  by  their  location  in  time  marked  on  the  dye 
curve  recordings  (Figure  12).  Usually  four  upstream  and  four  downstream 


O  w  w  UJ 


TIME.  SEC  ARBITRARY  ZERO 

Figure  12.  Dose  IV  dye  curve 


samples  were  analyzed  for  krypton  and  tritium  content.  If  these  two 
sets  of  samples  were  centered  at  their  respective  peaks,  the  average 
krypton-to-tr itium  ratio  of  each  set  and  the  lapse  time  from  peak  to 
peak  were  used  to  calculate  K^r  (dose  III).  If  a  set  of  samples  was 
centered  off  the  peak,  then  only  those  samples  centered  around  or  near 
the  peak  were  used  to  compute  K^r  .  In  this  case  the  temporal  center 
of  the  samples  used  to  compute  K^r  was  used  to  determine  travel  time 
(dose  II). 

53.  Samples  were  then  prepared  for  analysis  using  the  method 
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described  by  Cohen,  et  al.  The  krypton-85  and  tritium  concentrations 
were  determined  simultaneously  in  a  Packard  Tri-Carb  liquid  scintilla¬ 
tion  counter.  Three  replicates  of  each  sample  were  prepared  and  cycled 
through  three  counting  sequences  to  reduce  the  effect  of  analytical 
error. 


DO  Uptake  Tests 

54.  After  the  tracer  tests  had  been  completed,  tests  for  DO  up¬ 
take  using  the  disturbed  equilibrium  technique  were  conducted  in  the  hy¬ 
draulic  model.  Under  similar  flow  and  design  conditions  as  previously 
described,  a  sodium  sulfite  and  catalyst  solution  was  injected  into  the 
flow  upstream  of  the  model.  The  chemicals  combined  with  free  oxygen, 
effectively  creating  an  oxygen  deficit  in  the  water.  The  flow  was  re¬ 
aerated  as  it  passed  through  the  model  thus  decreasing  the  oxygen 
deficit . 

55.  DO  concentrations  were  measured  at  the  same  sampling  points 
as  in  the  radioactive  tracer  studies.  A  standard  YSI  DO  probe  was  used 
to  determine  the  DO  content.  In  addition,  a  modified  Winkler  titration 
technique  was  employed  to  spot-check  measurements. 
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PART  V:  RESULTS 


56.  Table  3  presents  krypton-to-tritium  ratios  for  the  four 

tracer  tests.  Doses  I  and  II  were  the  tests  with  the  flip  lip;  doses 
III  and  IV  were  the  tests  without  the  flip  lip.  The  ratios  whose  sam¬ 
ples  were  centered  around  or  closest  to  the  dye  peaks  were  used  in  the 
analysis  for  .  The  time  of  flow  was  determined  by  the  location  of 

these  samples  on  the  dye  curves.  Plates  1-4  show  the  dye  curves  and 
temporal  locations  at  samples.  The  time  of  flow  for  each  dose,  which 
was  used  in  computing  the  exchange  coefficient  ,  is  marked  "best" 

on  these  plates.  Plates  5-8  show  krypton-to-tritium  ratios  versus  time 
on  a  semilog  coordinate  system.  The  slope  of  the  line  is  the  krypton 
gas  exchange  coefficient. 

57.  Table  4  shows  the  samples  whose  krypton-to-tritium  ratios 

were  used  to  calculate  K,  for  each  dose.  The  travel  times  associ- 

kr 

ated  with  these  "best"  samples  are  also  presented.  By  applying  Equa- 

T 

tion  19  the  exchange  coefficient  for  oxygen  Kqx  at  the  ambient  water 
temperature  T(°C)  was  computed.  Applying  a  correction  for  temperature, 
the  reaeration  rate  coefficient  for  oxygen  at  20°C  was  determined  by  the 
following  equation: 


K20  =  KT  1.022(20-T) 
OX  ox 


(22) 


The  krypton  gas  fraction  remaining  in  the  water  at  the  downstream  sta¬ 
tion  was  computed  by 


where 


F.  R. 


r 


downstream 


r 


upstream 


F.  R.  =  fraction  remaining 


r 


upstream’ 


r,  =  mean  krypton-to-tritium  ratio  for  "best" 

downstream 

samples  at  the  upstream  and  downstream 
station 
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Dose  I 


F 

1.3415 

G 

0.8326 

G 

Not  available 

H 

0.8464 

H 

1.3563 

I 

0.8413 

I 

1.3438 

J 

0.8005 

Dose  II 

D 

1.2577 

B 

0.7854 

E 

1.2554 

C 

0.7486 

F 

1.2662 

D 

0.7338 

G 

1.3064 

E 

0.7314 

Dose  III 

F 

1.3694 

I 

0.9364 

G 

1.3811 

J 

0.9189 

H 

1.3684 

K 

0.9471 

I 

1.3370 

L 

0.8966 

Dose  IV 

D  1.2585  I  0.9568 
E  1.2650  J  0.9688 
F  1.2712  K  0.9158 
G  1.2765  L  0.9111 


I 


Table  4 

Samples  Used  in  K^r  Analysis  and  Travel  Time 


Dose 

Samples 

Upstream 

Included 

Downstream 

Travel  time 

sec 

I 

F,  H,  I 

I,  J 

22.7 

II 

F,  G 

B,  C,  D,  E 

24.5 

III 

F,  G,  H,  I 

I,  J,  K,  L 

28.5 

IV 

E,  F,  G 

J,  K,  L 

28.5 

For  example,  the  gas  fraction  remaining  at  the  downstream  station  for 
dose  I  was 


0.8209 

1.3472 


0.6093 


The  gas  fraction  lost  to  the  atmosphere  was 

1.000  -  0.6093  =  0.3907 

Similar  calculations  were  made  for  the  remaining  doses. 

58.  Table  5  summarizes  the  tracer  tests.  The  flip  lip  increased 

20°C 

the  gas  transfer  by  37.2  percent  and  increased  Kqx  by  80 . o  percent. 
Travel  time  decreased  by  17.2  percent. 

59.  During  dose  III,  air  bubbles  which  were  entrained  in  the 
water  were  being  pumped  through  the  sampling  system.  To  avoid  possible 
loss  of  tracer  gas  in  the  samples  due  to  these  air  bubbles,  the  pump  was 
positioned  slightly  lower  in  the  water  prior  to  dose  IV.  Repositioning 
the  pump  may  have  caused  the  difference  in  gas  loss  for  doses  III  and 
IV.  This  difference  was  larger  than  expected. 

60.  Figure  13  is  a  comparison  of  DO  predicted  with  the  and 

t  obtained  from  the  tracer  study  and  the  DO  observed  in  the  disturbed 
equilibrium  (DE)  tests.  The  small  variations  about  the  perfect  agree¬ 
ment  line  were  probably  due  to  incomplete  chemical  reactions  in  the  DE 
tests  or  slight  differences  in  the  mixing  conditions  in  the  model  for 
the  DE  tests  and  tracer  tests. 
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Table  5 
Test  Summary 


Design 

Dose 

Travel 

Time 

sec 

Fraction 

of 

Gas  Lost 

Water 
Tempera - 
ature 
°C 

„20°C  . 

K  per  hr 

ox 

With  flip  lip 

I 

22.7 

0.J9G7 

29.5 

77.0 

With  flip  lip 

II 

24.5 

0.4171 

29.5 

77.7 

Mean 

23.6 

0.4039 

77.3 

Without  flip  lip 

III 

28.5 

0.3220 

30.5 

47.1 

Without  flip  lip 

IV 

28.5 

0.2667 

29.5 

38.4 

Mean 

28.5 

0.2944 

42.8 
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PART  VI:  CONCLUSIONS 


61.  The  tests  showed  that  the  krypton-85/tritium  tracer  technique 
was  applicable  to  a  model  of  a  hydraulic  structure  for  determining  the 
gas  transfer  characteristics  of  a  particular  model  design.  The  tech¬ 
nique  provided  excellent  reproducibility  as  indicated  by  the  small  vari¬ 
ation  in  results  from  replicate  tests.  The  technique  was  demonstrated 
to  be  sufficiently  sensitive  to  be  used  as  a  tool  for  a  qualitative  com¬ 
parison  of  the  gas  transfer  of  the  different  designs  tested.  However, 
the  technique  may  not  be  sensitive  to  more  subtle  differences.  The  hy¬ 
draulic  model  which  was  tested  had  very  turbulent  hydraulic  conditions 
and  large  discharge  compared  to  other  hydraulic  models.  Thus,  there  may 
exist  a  threshold  below  which  sufficient  gas  transfer  does  not  occur  to 
permit  detection. 

62.  A  comparison  of  DO  observed  in  the  sodium  sulfite  tests  and 
predicted  DO  based  on  the  tracer  tests  was  quite  favorable  (Figure  13). 
Thus,  either  technique  would  provide  adequate  information  on  gas  trans¬ 
fer  for  this  particular  model  and  these  design  changes.  The  sensitivity 
threshold  for  the  DE  technique  may  be  different  from  the  tracer  tech¬ 
nique  threshold  because  of  the  different  measurement  and  analysis  tech¬ 
niques  employed.  Further  testing  is  necessary  to  evaluate  the  tech¬ 
niques  in  greater  depth. 

63.  In  a  recirculating  water  supply  system  such  as  that  used  at 
the  model  facility  described,  contamination  of  water  entering  the  model 
with  water  containing  chemicals  is  possible.  It  is  therefore  imperative 
that  the  circulation  time  be  much  greater  than  the  time  required  for 
data  collection.  The  sodium  sulfite  technique  is  not  feasible  if  the 
model  is  large  such  that  the  discharge  is  greater  than  0.25  cms,  since 
large  quantities  of  sodium  sulfite  solution  must  be  prepared  for  the 
high  injection  rate  necessary  for  this  large  flow. 
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PART  VII:  RECOMMENDATIONS 


64.  It  is  recommended  that  this  type  of  model  study  be  conducted 
for  projects  (proposed  or  existing)  where  modifications  are  contemplated 
to  improve  gas  transfer.  A  qualitative  evaluation  of  different  struc¬ 
tural  alternatives  may  result  in  a  more  effective  or  economical  design. 
It  is  probable  that  prototype  response  is  directly  related  to  changes  in 
gas  transfer  characteristics  observed  in  a  model.  However,  relation¬ 
ships  for  quantitative  transfer  of  gas  exchange  characteristics  from 
model  to  prototype  have  not  been  developed  yet.  Increased  gas  exchange 
in  a  model  directly  indicates  greater  oxygen  transfer  to  low  DO  water  in 
a  prototype.  Increased  gas  exchange  in  the  model  should  indicate  de¬ 
creased  supersaturation  if  the  prototype  causes  supersaturation.  It  is 
also  recommended  that  similar  tests  be  conducted  in  hydraulic  models 
that  range  in  size  and  scale  to  determine  the  gas  transfer  occurring  in 
other  hydraulic  models  and  determine  the  effects  of  other  models  on  the 
applicability  of  the  technique.  Other  types  of  simple  geometry  should 
be  tested  to  provide  a  basis  for  scaling  criteria  and  to  evaluate  detec¬ 
tion  thresholds  below  which  the  technique  is  not  adequate. 

65.  Research  to  determine  scaling  relationships  should  be  con¬ 
ducted  using  physical  models  of  simple  turbulent  flow  systems.  A  param¬ 
eterization  study  to  isolate  and  evaluate  the  significant  parameters  af¬ 
fecting  the  scaling  of  gas  transfer  in  models  should  be  conducted. 
Theory,  applicable  to  scaling  gas  transfer,  should  be  developed  concur¬ 
rently  with  the  parameterization  study.  Thus  an  interface  is  provided 
for  evaluating  experimental  and  theoretical  results.  This  approach  is 
being  pursued  at  the  Waterways  Experiment  Station  in  research  funded  by 
the  Environmental  and  Water  Quality  Operational  Studies.  It  is  hoped 
that  this  will  result  in  a  method  of  estimating  scale  relationships  be¬ 
tween  the  gas  transfer  characteristics  of  a  model  and  prototype  and  the 
prediction  of  oxygen  uptake  in  the  prototype. 
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PLATE  8 


DOSE  IV 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Wilhelms,  Steven  C 

Tracer  measurement  of  reaeration:  Application  to  hydraulic 
models;  hydraulic  model  investigation  /  by  Steven  C.  Wilhelms. 
(Hydraulics  Laboratory.  U.S.  Army  Engineer  Waterways  Experiment 
Station)  ;  prepared  for  Office,  Chief  of  Engineers,  U.S.  Army, 
under  CWIS  310l*2,  (EWQOS  Work  Unit  3l6ol*  (IIIA.2)  —  Vicksburg, 
Miss.  :  U.S.  Army  Engineer  Waterways  Experiment  Station; 
Springfield,  Va.  :  available  from  NTIS,  1981. 

p.  ,  t1*]  leaves  of  plates  :  ill.  ;  27  cm.  (Technical 
report  /  U.S.  Army  Engineer  Waterways  Experiment  Station  ; 

E-80-5 ) 

Cover  title. 

"December  1980." 

References:  p.  bl— 3 . 

1.  Hydraulic  models.  2.  Radioactive  tracers.  3.  Reaeration. 

1* .  Tracers.  I.  United  States.  Army.  Corps  of  Engineers.  Office 
of  the  Chief  of  Engineers.  II.  United  States.  Army  Engineer 
Waterways  Experiment  Station.  Hydraulics  Laboratory.  III.  Title. 
IV.  Series:  Technical  report  (United  States.  Army  Engineer 
Waterways  Experiment  Station)  ;  E-80-5. 

TA7.W31*  no.  E-80-5 


